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A renormalization group theory for a system consisting of coupled superconducting layers as a 
model for typical high-temperature superconducters is developed. In a first step the electromagnetic 
interaction over infinitely many layers is taken into account, but the Josephson coupling is neglected. 
In this case the corrections to two-dimensional behavior due to the presence of the other layers are 
very small. Next, renormalization group equations for a layered system with very strong Josephson 
coupling are derived, taking into account only the smallest possible Josephson vortex loops. The 
applicability of these two limiting cases to typical high-temperature superconductors is discussed. 
Finally, it is argued that the original renormalization group approach by Kosterlitz is not applicable 
\ to a layered system with intermediate Josephson coupling. 
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I. INTRODUCTION 



Since the highly anisotropic high-temperature super- 
conductors (HTSC's) consist of weakly coupled super- 
conducting layers one expects two-dimensional effects 
to be important. Quasi- two-dimensional superconduct- 
' ing films near the critical temperature and in the ab- 
sence of an external magnetic field are well described 
■ by the Berezinskii-Kosterlitz-Thouless (BKT) theory.El 
The BKT theory states that the unbinding of sponta- 
neously created vortex-antivortex pairs is responsible for 
O \ the phase transition of these films. A rigorous renormal- 
O ■ ization group .formulation of the BKT theory was given 
^ ' by Kosterlitz. □ The question arises of whether the BKT 
theory can be extended to describe layered superconduc- 
^ ' tors. 

In this paper we first consider a system of supercon- 
ducting layers with exclusively electromagnetic coupling 
and than extend our approach to include the Josephson 
coupling between the layers. The main consequence of 
the Josephson coupling is the appearance of Josephson 
vortices (JV's), i.e. vortex strings that reside between the 
superconducting layers (for a review see Blatter et al.u) . 

The layered system without Josephsoii, coupling has 
been treated by Scheidl and Hackenbroicba within a self- 
consistent linear response theory. The authors describe 
the renormalization of the interaction between vortices 
residing in the same layer, in neighboring layers, etc., by 
a single equation. It is shown below, however, that the 
screening cannot be described in such a simple way. 

Sewegal authors-have considered the Josephson coupled 
case-Era Horovitza has studied the phase transitions due 
to a) the unbinding of vortex-antivortex pairs with van- 
ishing Josephson coupling and b) JV loops that reside en- 
tirely between the layers within a sine-Gordon renormal- 
ization group approach. The competition between these 
two mechanisms is then discussed in a partly heuristic 
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manner. The role of JV loops that cross the layers does 
not become clear. 

FischerQ has assumed that the renormalization of the 
linear term in the interaction of pancake vortices, which 
is due to the presence of JV's, and the logarithmic, elec- 
tromagnetic term can be described by the same screening 
factor. This is not the case. 

PiersonQ has used a renormalization group approach 
similar to Ref. |2| to treat the Josephson coupled system. 
However, the partition function given by Pierson does 
not correctly describe a layered superconductor because 
it contains a linear term in the vortex interaction for any 
two pancake vortices in the same layer or in neighboring 
layers, whereas such a term should only be present for 
two pancakes that are actually connected by a JV. 

Friesenta has considered the further simplified model of 
one superconducting layer between two thick supercon- 
ducting slabs, using an expression for the vortex inter- 
action which is correct for small separations only. This 
approximation is not sufficient since in the renormaliza- 
tion procedure the interaction between widely separated 
vortices must be taken into account. 

A different approach to the same jproblem has been 
used by Chattopadhyay and Shenoy.EI They start from 
the anisotropic three-dimensional (3D) XY model on a 
cubic lattice. There are two reasons, however, why we 
prefer to use a different model: 1. The HTSC's are es- 
sentially discrete in one direction, but continuous in the 
other directions, whereas the XY model of Ref. ^ is dis- 
crete in all directions. Universality ensures that both 
systems behave similarly within the 3D critical region, 
but the present paper is concerned with the 2D behav- 
ior outside that region, where universality does not hold. 
2. The XY model on a cubic lattice contains three in- 
dependent parameters (two spin couplings and a lattice 
constant), whereas the model considered here contains 
five parameters in the case of non-vanishing Josephson 
coupling. Experiments suggest that there are no uni- 
versal relations between these parameters, i.e. they are 
truly independent. Thus the 3D XY model appears to 
be insufficient to describe the HTSC's. Furthermore, like 
PiersonD and FriesenEl the authors make the approxima- 
tion of a linear term in the interaction between any two 
pancake vortices, irrespective of the positions of the JV's. 

The present paper is organized as follows: In Sec. |n] a 
rigorous renormalization group theory for the exclusively 
electromagnetically coupled system is developed and the 
quantitative predictions of this theory are discussed for 
a typical HTSC. This theory then serves as a reference 



frame for the further discussion. In Sec. |HJ the con- 
sequences of Josephson coupling are investigated. The 
paper concludes with a brief summary. 



II. THE ELECTROMAGNETIC 
APPROXIMATION 
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A. Derivation of the renormalization group 
equations 



Let us now turn to the layered system with negligible 
Josephson coupling. This approximation may be appli- 
cable to artificially grown HTSC superlattices since these 
structures are characterized by a large layer separation 
compared to the coherence length perpendicular to the 
layers. 

The unrenormalized interaction between two vortices 
separated by a distance r within the same layer is given 
byl3 



Uo(r) = T^§^h, r -, (2.1) 



r 



where the upper (lower) sign is for vortices of the same 
(opposite) vorticity, 0o = hc/2e is the flux quantum, s 
is the layer separation, A a {, is the magnetic penetration 
depth for fields applied perpendicular to the layers, and 
r is the minimum vortex separation. The interaction 
between vortices that are n ^ layers apart isO 

Note the opposite sign and the reduction by the small fac- 
tor s/2A a h compared to the interaction within the same 
layer. Note further that the interaction of vortices that 
are, say, 2 or 3 layers apart is of the same order of mag- 
nitude as the interaction for n = 1. Thus it would not be 
justified to include only the interaction between vortices 
in neighboring layers and neglect terms with \n\ > 1. 
Defining the vortex "charge" 



9=± V^ (2 ' 3) 



and coupling parameters 



— 1 for n = 

3 f \n\s\ (2.4) 

— — exp — - — I for n f= 

2A Q b \ X a b J 

we can write down the grand canonical partition function, 

'■Eipt/E/^-E/, 



N Hi 1 U2N 

x exp 



- 2^ miO-m-nj In I • (2.5) 

Here, N is the total number of vortex-antivortex pairs 
in the system, z = exp(/3/x) is the vortex fugacity, 
f3 = 1/fegT is the inverse temperature, n, and are the 
number of the layer and the position within the layer of 
vortex i, respectively, and the range of integration Di is 
the whole of the layer rij except discs of radius r centered 
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at the vortices j < i. Since the external magnetic field 
vanishes the sums over the layer indices Hi are subject to 
the constraint that the total vorticity in every layer be 
zero. 

Following the program of the renormalization group 
approach,B we now proceed to integrate out the smallest 
vortex-antivortex pairs (of size between r and t + dr). 
The central concept is to rewrite the resulting expres- 
sion in such a way that it has the same functional form 
as the original partition function, but with renormalized 
parameters. In our case, these parameters are the fu- 
gacity z and the coupling parameters a n . As more and 
more pairs are integrated out, the parameters are renor- 
malized to incorporate the effect of these pairs on the 
system. The renormalization is described by differential 
equations for the parameters as functions of the smallest 
pair size r. The properties of the macroscopic system 
are governed by the limiting behavior of these recursion 
relations for large t. 

In the case of a single layercl only neutral pairs are in- 
tegrated out so that the total vorticity remains zero. A 
change in the total vorticity would result in an infinite 
contribution to the energy and could not be compensated 
for by suitable renormalization of parameters. Similarly, 
in the present case the vorticity in every layer must re- 
main zero because removal of any configuration of vor- 
tices that does not consist of pairs in the layers would 
make the energy of the remaining system infinite. (This 
observation is related to the fact that the interaction is 
screened by vortices which move freely within the planes, 
but not perpendicular to the planes.) Thus only pairs 
which reside in the same layer can be integrated out. 
This limitation of the present approach is expected to 
have unphysical consequences, especially concerning the 
correlations between pancakes in different layers, which 
are essential within the (narrow) 3D critical region. 

The integration procedure is similar to the case of a single layer .□ It may be written as the following prescription: 



E / <*V"E / <?r 2N = J2 / /rvE / , 

m n 2N J D 2N ni JD[ „ 2JV JD' 2N 

^EE E L E L <*w-E /„. d ViE L ^ 



ijkj ni u D i nt-i D i-i n i+1 D i+i rij-i D j-i n j+i D i+i 



j+1 



2« nj Jd j m Jt< - 



x } d r 2N > \_d rj ) d n6 (2.6) 

r< \ Ti — Tj <T+dT 



Here, the D[ are the same as the Di above, but with r replaced by r + dr, and Dj is the whole layer except discs 
of radius r centered at the vortices k ^ The sum over can be performed trivially. The right-hand side of 



Eq. (2.6) consists of two summands: The first one is similar to the whole left-hand side, but with minimum pair size 



r + dr instead of r. The second one gives the approximate correction. The last integral in Eq. (2.6) integrates over 
the separation vector of a small pair. The integral over rj and the sum over rij take the pair over the whole layer and 
all the layers, respectively. The sum 1/2^^ - selects every possible pair just once. 



We now proceed to apply the prescription (2.6) to the partition function (2.5). The mathematical procedure is 
similar to the one of Ref. ^, taking into account the additional sums over layer indices. The result is 



exp 



27r( — I rdrMF 



d r 2 N 



4 



x exp 



2 ^ 



dr \ 

tm-nj + 2ir 2 z 2 —f3q 2 a ni _ nj J 



In 



r« - r , 



(2.7) 



where <S„ = ^ m a m _„a m and -F and M are the size and the number of the layers, respectively. Up to this point, r 
has been replaced by r + <ir only in the ranges of integration D\. However, we have to rescale r everywhere to be 
consistent. To order dr we eventually obtain 



exp 



7T (4-) TdrMF 
d 2 n 



i 



^ AH 2 \{T + dr) 2 



2N r 



E 



d r 2 A?exp 



2 SJ 



and 



2n 2 z 2 — /3q 2 a n . 

T 



dz ( _ ft 2 



da,, 



2ir 2 z 2 Pq 2 a n , 



where a logarithmic length scale I — Itlt/tq is introduced. 
Here, ro is the unrenormalized minimum pair separation. 
These equations may be rewritten as 



dz 2 

— = z 2 (4 + f3q 2 a ), 
dl 

d a n „ 2 2 a 1 

-— = 2tt z [jq y a m _ K a m . 



dr 

2Tr 2 z 2 —Pq 2 a nz . 

T 



-i 2N 



Note that no length parallel to the z axis is rescaled, 
in contrast to the length scale in the planes, r. Indeed, 



the partition function (2.5) does not contain any length 



scale in the z direction. We can interpret the layer index 
as an internal, discrete degree of freedom. The system 
we get in this way is equivalent to the layered model 
since both have the same partition function, but it is two- 
dimensional. Thus, rescaling of lengths in the z direction 
is meaningless for both systems. 

If we comp are Eq. ( |2.S| ) with the original parti- 
tion function (|2.5|), we find that both have indeed the 
same functional form. Dropping the irrelevant factor 
exp(27r(z/r 2 ) 2 r(irMF), the partition function takes the 
form of the original partition function if we set 



(2-9) 



(2-10) 



The same information can be expressed as a system of 
coupled differential equations: 



(2.H) 
(2.12) 



(2-13) 
(2.14) 



In- 



dr 



(2. 



This infinite set of equations replaces the Kosterlitz re- 
cursion relations in the case of an electromagnetically 
coupled layered system. It is subject to the boundary 
condition that z 2 and a n take on their bare values for 
/ = 0. Note that by letting ao = — 1 and a n — for 
n ^ we regain the Kosterlitz recursion relations. 
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B. Weak electromagnetic coupling 



The macroscopic properties of the layered supercon- 
ductor are contro lled by the behavior of the solution of 
Eqs. ( 2.13 ) and ( 2.14 ) for large length scales I. These 
equations can be simplified by means of an expansion 
for small electromagnetic interaction between vortices in 
different layers. 

To this end we define the quantity 



A 2 



On 



(2.15) 



which is a measure for the coupling between different lay- 
ers as compared with the coupling within the same layer. 
From the definition (|2.4|) it follows that in the unrenor- 
malized case A 2 = s/4X ab < 1 for all HTSC's. We will 
see below that this inequality holds for the renormalized 
quantities also. 

If we further define the usual stiffness constant 



K(l) = 



ao(l)Pq 2 
2tt 



we find, to linear order in A 2 , 
2z 2 (2-ttK), 
-Att 3 z 2 K 2 (1 
= -8tt 3 z 2 KA 2 . 



dz^ 

If 

dK 
~dT 
dA 2 



A 2 



dl 

Formal integration of the last equation yields 

cl 



A 2 {I) = A„exp 



8tt 3 / dl' z 2 (l')K{l') 



(2.16) 

(2.17) 
(2.18) 
(2.19) 

(2.20) 



where A^ 



A 2 (l = Q). 



Since A^ ^ 1 we may expand the square of the fugacity 



and the stiffness constant to linear order in Ak 



0' 



z 
K 



zl + Az 2 Al, 



K 



AKAl- 



(2.21) 
(2.22) 



To the same order, z 2 and K in Eq. (2.20) can be re placed 
by z 2 and K . Inserting the expansions into Eqs. ( 2.17 ) 
and (2.18), we obtain a set of four coupled equations, 



dl 
dAz 2 
dl 

dKo 
dl 
dAK 

dl 



= 2z 2 (2~ttK q ), 
= 2(2-TrK Q )Az 2 
= -4tt 3 z 2 K 2 , 



2nz 2 AK, 



(2.23) 
(2.24) 
(2.25) 



-4n 3 K 2 Az 2 - 8ir 3 z%K AK - Att 3 z^K 2 



x exp 



8tt 3 / dl' z 2 {l')K (l') 



(2.26) 
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These equations are subject to the boundary conditions 
4(0) = z 2 {0) = exp(2/3/z), Az 2 (0) = 0, K (0) = K(Q) = 
-an(O)0Q 2 /2TT, and A£[0) = 0. With the help of 



Eq. (|2.25D the integral (|2.20|) is found to be 

A2 v =A2 °{w)) 2 - (2 - 27) 

Thus we obtain a purely differential equation for AK: 
"' AA ' - -^K 2 Az 2 - %^z 2 KqAK - A^z 2 K ° 



(2.28) 



Equations fl2.23|) to fl2.25| ) and Eq. ( p.28|) describe the 



layered system for small s/4A a (,. Note that the num- 
ber of equations has be en re duced from infinity to four. 
Note further that Eqs. ( |2.2% and ( |2.25| ) are the original 



Kosterlitz recursion relations p which describe the uncou- 
pled layers and can be solved by themselves. 

The solution of the Kosterlitz recursion relations is well 
known. For T smaller than a critical temperature T c 
the fugacity zq converges exponentially to zero for large 
I, while the stiffness goes to a finite value Kq > 2/vr. 
Thus the vortices are bound in small vortex-antivortex 
pairs. For T > T c the fugacity diverges exponentially and 
Kq goes to zero for I — ► oo. Here, very many unbound 
vortices exist. At T c the stiffness constant Kq jumps from 
2/-7T to zero; this is the famous "universal jump" . 

Let us now take the electromagnetic coupling between 
the layers into account. The asy mptot ic behavior of 



the quantity A 2 is given by Eq. ( [2.27 ), whereas the 



asymptotic forms of Az 2 and AK can be read off from 



Eqs. (2.24) and (2.2S) in connection with the known 
forms of z'q and Kq. 

We first consider the case T > T c . The quantity 
Az 2 diverges exponentially for large Thus the fu- 
gacity z — \/ Zq+ Az 2 Aq also diverges exponentially, 
whereas AK, K, and the coupling parameter ao van- 
ish. Thus the fugacity and the stiffness constant qualita- 
tively behave as in the case of a single layer. Note that 
d/dl(Az 2 j ' Zq) = —2nAK so that Az 2 / z 2 app roaches a 



constant for I — > oo and the expansion (2.21) remains 
valid for arbitrarily large I. From Eq. ( p. 27 ) it follows 



that A 2 goes to zero in this regime. Since A z is a measure 
of the electromagnetic coupling between the layers this 
means that the renormalized interaction between vortices 
in different layers vanishes faster than the interaction 
within the same layer. In this sense the layers decou- 
ple and the system becomes two-dimensional for T > T c . 

For T < T c , on the other hand, the quantity Az 2 and, 
thus, the fugacity z converge towards zero. The correc- 
tion to the stiffness, AK, approaches a finite and negative 
value. Therefore the coupling parameter ao of vortices 
within the same layer is reduced by the presence of other 
layers, the correction being proportional to s/X a b- As 
long as this correction is small, the pairs are still bound, 
though less tightly. The quantitative significance of this 
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correction is discussed below. The quantity A 2 is reduced 
by the renormalization, but remains finite. Thus the in- 
teraction between vortices in different layers remains fi- 
nite, and the system is three-dimensional for T < T c . 

Whereas the above results are obtained by analyti- 
cal study of the asymptotic behavior of the recursion 
relations, numerical integration is necessary to produce 
trajectories in the AK-Az 2 plane. Trajectories for 
z(0) — 0.03 and several starting values K (0) are shown 
in Fig. |l|. The value z(G) = exp(/3u) = 0.03 seems 
reasonable for typical HTSC's (Ref. || gives the value 
fi = —E c = — 0.72g 2 and [3q 2 must be larger than 4). 

We now turn to the question of the absolute size of the 
correction AK to the stiffness constant below the critical 



temperature. Numerical studies of Eqs. ( p.24| ) and ( 2.2S ) 
show that 

, N 27r 3 z 2 (0)K n 2 (oo) 

is a very good approximation for small z(0) and small 2 — 
TrK a (po). The temperature dependence of K a is known 

tob±3 



K (oo)Si- [I + VWJ^— , (2.30) 




where B is a nonuniversal constant. Thus 

Aif(oo) S -27rz 2 (0)\/— - (T c - T)~ 1//2 (2.31) 
V B 

to leading order in T c — T. 

If the critical temperature is approached from below, 
the correction to the stiffness constant thus diverges with 
an exponent of —1/2. Very near to the transition tem- 
perature of a single layer, T c , this divergence causes the 
stiffness constant K(oo) to become negative, which is 
physically impossible. In this region higher order terms 
neglected here should render K > 0. We expect that K 
vanishes and the pairs start to break up at a temperature 
T C 3D < T c . T 3D should be larger than the temperature T* 
at which K vanishes in the linear approximation (which 
thus defin itely fails at T*). T* can be estimated from 



Eq. ( [2.31D . For Bi-2212 the parameters are z(0) » 0.03, 
s/4A a6 w 0.0019, and V2B » 2.9lH We thus obtain 
T*/T c = 1 - 1.3 • 10" 10 . Obviously, T* is indistinguish- 
able from T c , the linear approximation is always valid 
in practice, and the transition temperature T 3D is not 
significantly shifted by the electromagnetic coupling. 

The foregoing discussion indicates that the effect of the 
other layers is very small. Indeed, even if we could mea- 
sure a correction of AK/Kq = 1% in Bi-2212, we would 
still need to resolve a temperature range of the order of 
10~ 6 T C to see it, which is impossible. In practice the ef- 
fect of the electromagnetic coupling will be overruled by 
the Josephson coupling, see below. 

By using HTSC superlattices it may be possible to in- 
crease the ratio s/X a b- An additional advantage of su- 
perlattices is that the electromagnetic approximation is 
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more appropriate in this case. However, even then the 
effect will be very small. More promising systems are 
superlattices fabricated from conventional superconduc- 
tors. 



III. LAYERED SUPERCONDUCTORS WITH 
JOSEPHSON COUPLING 

A. General consequences of Josephson coupling 

Since Josephson tunneling between the layers has been 
found experimentallyEi in Bi-2212 and other HTSC's 
we are faced with the question how to incorporate the 
Josephson coupling into the renormalization group the- 
ory presented above. 

The Josephson coupling between the layers leads to 
the appearance of JV's. Because of the smoothness of 
the phase of the order parameter in the layers and the 
conservation of magnetic flux every two-dimensional vor- 
tex (i.e. pancake vortex) in a layer must be connected to 
two JV strings and vice versa. Since the external mag- 
netic field vanishes and surface effects are neglected the 
JV's form vortex loops. On these loops pancake vortices 
sit, where ever the loops penetrate a layer. It is clear 
from this picture that a JV connects either a vortex and 
an antivortex within the same layer or two vortices of the 
same vorticity in neighboring layers. 

The program of the renormalization group approach to 
this particular system is to derive an expression for the 
energy of any possible configuration of JV's, add this ex- 
pression to the Hamiltonian, insert the new Hamiltonian 
into the partition function, integrate out small vortex- 
antivortex pairs, and derive a set of recursion relations. 
To obtain a tractable expression for the energy of JV's 
we introduce several approximations. 

The energy of JV strings is_aporoximately proportional 
to their length L for L > AjJijliZj where Aj is the Joseph- 
son length. If two pancakes are separated by less than Aj, 
no full JV develops and the estergy due to the Josephson 
coupling is not simply linear However, for small sepa- 
rations the electromagnetic contribution to the interac- 
tion dominates and the form for the Josephson term is 
unimportant. Furthermore, as the interaction is renor- 
malized, small pancake pairs are integrated out and only 
large pairs remain, for which the linear approximation 
is correct. Thus we assume that the energy of a JV of 
length L is given by 

Uj = kL. (3.1) 

It is further assumed that the JV's form straight lines 
between the pancakes. Thus we may take the JV's into 
account by adding a term of the form ([O]) to the Hamil- 
tonian for any two vortices that are connected by a JV. 
The Hamiltonian takes the form 
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H = \ E '/''(/°- " ln Vl ^ + E K lv- r pL 

' (P:P')ec./ 

(3.2) 

where Cj denotes a given configuration of JV's connect- 
ing the pancakes i = 1, . . . , 2N and (p,p') is a single JV 
string characterized by the pancakes p and p' at its ends. 



The constant k in Eq. (3.2) can be evaluated for small 
s/^ab, the unrenormalized value beingjlj k(1 = Q) = 
<Z 2 A a b/sA c , where A c is the penetration depth for mag- 
netic fields applied parallel to the layers. Note that 
Kb/sX c « A -1 for Bi-2212. We can now see that 
the electromagnetic approximation is not sufficient for 
Bi-2212: The average distance between neighboring pairs 
and the separation on which the linear term in the inter- 
action becomes important are of the same order of mag- 
nitude (1000 A). 

Since the partition function is a sum over all possible configurations, we must, in addition to the summations and 



integrations in Eq. (2.5), sum over all configurations of JV's for given pancake positions. Thus the grand canonical 

(3.3) 



partition function may be written as 

2JV 



Z = E j& (£) " E / • • • E / *™ £ eM-W) 



with TL given by Eq. ( |3.2| ). The last sum in Eq. (3.3) is over all con figurations Cj of JV's. 

We are now faced with the task to rewrite the prescription ( |2.6| ), which tells us how to integrate out small pairs, 
for the partition function ( |3.3| ). The new aspect here is that we must take all possibilities into account to insert the 
small pair into the JV loops. These possible insertions are of three topologically different types: 1. the small pair can 
form a vortex loop by itself, as shown in Fig. |](a). 2. the pair can be part of a larger vortex loop, but still be directly 
connected by one JV, see Fig. ||(b). 3. the pair is not directly connected by a JV, see Fig. ||(c). Sin ce th e present 
approach neglects some contributions to the interlayer vortex correlations as discussed in subsection 



Q A 



we may 

expect unphysical results if vortex loops across more than one layer are essential. On the other hand, the approach 
is consistent if only small loops of type 1 are important . 

If the sum over all possible insertions of the small pair into the configuration Cj of JV's between the other 
pancakes is denoted by 



E 



(3.4) 



the prescription for integrating out small pairs is given by 



E / ^-E / ^E=E / «*v-E / ^-E 

»U JDl n 2N ° 2N Cj »U JD 'l n 2N D 2N Cj 

ijij ni JD 'l n l - 1 jD i-l n i + 1 jD ', + l n J ^ 1 jD 'j-l n j+1 J D 'j+1 



2 r j+1 ■ ■ ■ 



E / ' /2 '-^v EE /_ 1,1 r j E / ,,2r ' E s ^ s qi ,- qj . ■ ■■!.■. ) 

n 2N jD 2N Cj n 3 J Dj m " , T<|r I -r J |<r+dr {l ^ C 'j 



Here, Cj denotes a configuration of JV's if the pancakes i and j are absent. Note that the sum (3.4) contains one 
summand for the first type of vortex loops and many summands for the second and the third type. Before the general 
case is discussed, we turn to a simple limiting case. 



B. Very strong Josephson coupling 

Let us now consider a layered superconductor with very strong Josephson coupling between the layers. In this case 
the energy of the JV strings is generally large compared to the other energy scales. Fluctuations of the JV's cost 
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much energy and the system prefers states with small total JV length for any given configuration of pancake vortices. 
Although states with perfectly aligned pancakes in all layers have the lowest Josephson energy, such states are not 
created by thermal fluctuations since their total energy is infinite. The elementary excitations are small vortex loops 
across one layer (see Fig. ^|(a)). Larger, more complicated loops will not form, except within a narr ow 3D region, 
since the energy of JV strings is large. Thus we consider only loops of the first type, and the sum (3.4) is reduced to 
one term. 

The limit of large k is not physically accessible since for a layered superconductors one always finds A c > X a b- 
Nevertheless it is considered here because it forms the opposite limiting case as compared with the system without 
any Josephson coupling and should allow one to estimate the maximum effect Josephson coupling can have. 

The additional energy resulting from the two JV's is 2Uj = 2k |r< — r ^ | = 2kt. This is a constant term in the 
Hamiltonian. Therefore, a constant factor e ~ 2 ^ rtT appears in the integrand. 

For the same reasons as in the exclusively electromagnetic case only pairs within the same layer can be integrated 
out. The calculations are similar. The constant factor mentioned above is just dragged through. After the integration 
and the replacement of r by r + dr, we obtain the partition function 



Z = exp 



>H0 



2vr(-^) TdTMFe- 20KT 

T 



y- 

^ NY- 



8 , \dr 



2N 



N\ 2 \(T + drr¥ 



2.Y 



x exp 



2 ^ 



On* —n 



dr 

2ir 2 z 2 — 3q 2 a r , 

T 



Qiqj In 



t + dr 



2(5n \ r p'- r p\ 



( P ,p')eCj 



(3.6) 



Again, there is no rescaling of lengths in the z direction, in contrast to Refs. ^ and ^| However, the arguments given 
in subsection II A remain valid in the presence of Josephson coupling and, therefore, the system is still essentially 
two-dimensional. 



Comparison with the original partition function (3.3) 
yields the recursion relations 



^- = z 2 {A + f3q 2 a ), 
da 



—— = 2tt z fiq exp(-2/3KT e ) 2_^a m _ n a m , 



dl 
dn 



(3.7) 
(3.8) 

(3.9) 



where again I = lnr/ro- We thus find that the linear 
coupling constant k is not renormalized due to screen- 
ing if we consider only loops of the first type. (Neither 
is there renormalization due to rescaling of lengths, as 
noted above.) These equations, especially Eq. fl3.9|), will 
break down very near to the phase transition due to the 
appearance of larger, three-dimensional vortex loops. 

The additional factor exp(— 2/?KToe') in the equation 
for a n approaches zero very rapidly. The change of the 
electromagnetic coupling parameters a n with the scale I 
and, therefore, the screening of the electromagnetic inter- 
action are strongly suppressed. In connection with the 
observation that the linear interaction due to JV's is not 
screened at all within this approximation, this result in- 
dicates that the low-temperature phase of bound pairs is 
stabilized and T c is higher than in the electromagnetic 
case. 



(3. 



To address this issue the recursion relations (3.7) and 
are expanded for small s/4A a h. The procedure is 



analogous to the one employed above and is not repeated 
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here. It is seen that the electromagnetic coupling between 
the layers brings about a very small correction to the 
renormalized electromagnetic coupling within the same 
layer. Thus it suffices to consider the vortex fugacity zq 
and the stiffness constant Kq to order zero. The renor- 
malization of these quantities is described by the equa- 
tions 

^ = 2z 2 (2-^ ), (3.10) 
^ = -4tt 3 z 2 X 2 exp(-2/3«T 0e '). (3.11) 

Note that f3q 2 = 2nK(l = 0). Numerical integration for 
2(0) = 0.03, k = 10~ 3 A~V, r = 21.5 A, and several 
values of K(0) yields the trajectories shown in Fig. 0. 



By studying the asymptotic behavior of Eqs. (3.10) and 
( |3.11 ) two temperature regimes can be identified. In this 



section the temperature which divides these two regimes 
is denoted by T c , whereas the original BKT transition 
temperature is Tbkt ■ As is shown in Fig. ^, T c is signif- 
icantly higher than Tbkt- For Bi-2212, experiments oa. 
ultra-thin films have found a transition temneraturecj 
Tbkt = 35 K and a mean-field temperaturelI3 T c q = 
86,8 K. From the above results, the transition temper- 
ature for bulk Bi-2212 shoukLbe T c = 37.4 K, but the 
experimental value is 84.7 K.Ej Taking into account that 
a more realistic description of the Josephson coupling 
(inclusion of larger loops) would probably yield an even 
lower T c , the main effect must be brought about by other 
mechanisms, e.g. change of the hole density in the CuC>2 
planes. 

For T < T c the fugacity goes to zero and the stiff- 
ness constant approaches a finite value Kq(oo) > 2/ir as 
I — > oo. The linear coupling due to JV's is not changed. 
The vortices are bound in small vortex-antivortex pairs. 

For T > T c the fugacity diverges exponentially, but 
the stiffness constant Kq still approaches a non- vanishing 
value i^o(oo) < 2/tt. Very many vortices exist, but they 
are still bound in pairs. This result is different from the 
exclusively electromagnetic case. Also, the electromag- 
netic interaction between vortices in different layers does 
not vanish and the system remains three-dimensional 
even above T c . 

Furthermore, numerical studies indicate that there is 
no special feature in K (cc) at T c (see Fig. ||). Nor do we 
expect the appearance of ohmic resistance at T c since the 
vortices are bound in pairs below as well as above this 
temperature. Thus the question arises of whether there 
really is a phase transition at T c . However, since the true 
transition is governed by effects not included here (larger, 
3D loops), this result should be interpreted with caution. 



C. Intermediate Josephson coupling — failure of the 
renormalization group approach 

Let us finally consider a layered superconductor with 
general Josephson coupling. In this case we must take 
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all the possible configurations of JV's into account (see 
Figs. ^|). This is especially necessary in the vicinity of T c . 
Here, we consider only one term that shows in what way 
the approach fails. 



It suffices to investigate one summand of the sum (3.4) 
of the second type (see Fig. ^(b)). Let ig (rj) be the 
position of the neighbor of vortex i (j) other than vortex 
j (i), see Fig. [|. The energy of these three JV's is given 
by k \ri — Tq| + kt + k \rj — rj|. Furthermore, to insert 
the new pair we must remove the J V from r l to r \ . Thus 
we have to consider an additional factor 

exp (-/3k \ri - r l Q \ - (3kt - 0k \rj - rj| + [3k \y{ - r^|) 

(3.12) 

in the integrand. The integral over the lowest order 
term in (3 of the factor containing the electromagnetic 
interaction — this term is simply unity — times the above 
factor is evaluated in the appendix. After integration 
over Vi and Yj we obtain the expression 

2-kt(It exp (^-(3kt + (3k \v[ - I q ((3kt) 

x|rj-r*| 2 X 2 (/3K|r{-r*|), (3.13) 

where Iq and K 2 are Bessel functions. Expansion for 
large separations r = |r{ — rg| yields 



2nTdTe-' }KT Io(p K T) ] J^-r^ 2 . (3.14) 

After performing the sums over the layer indices rii and 
rij and the vortex indices i and j we end up with a par- 
tition function based on a Hamiltonian that contains a 
term proportional to |ry — r p | 3 / 2 for any two pancakes p 
und p' that are connected by a JV. 

It can be shown that this term cannot be canceled 
by any other terms: All other terms in the integrand 
either contain factors (3q 2 or (3 2 q 4 (as opposed to the 
term discussed above) or have the same sign. 

But no term proportional to the vortex separation to 
the power 3/2 is present in the original Hamiltonian! 
Thus it is impossible to renormalize the parameters in 
the partition function in such a way as to regain the orig- 
inal partition function ( p.3j ). The program of the renor- 
malization group approach thus fails for the layered su- 
perconductor with (general) Josephson coupling, at least 
within the present model. We speculate that the reason 
for this failure lies in the one-dimensional nature of the 
JV strings as opposed to the point-like pancake vortices. 



IV. CONCLUSIONS 

The application of the cenormalization group formu- 
lation of the BKT theoryo to layered superconductors 
has been investigated in detail. This program has suc- 
cessfully been carried out for a layered system without 
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Josephson coupling. It has been found that the transi- 
tion temperature T c is not affected by the presence of 
other layers and that the correction to the stiffness con- 
stant K is significant only in a inaccessibly narrow tem- 
perature range below T c for typical HTSC's, but may be 
observable for layered structures made of conventional su- 
perconductors. The opposite limiting case of very strong 
Josephson coupling has also been investigated. A signifi- 
cant upward shift of the transition temperature has been 
found, but even using approximations valid only for un- 
physically strong Josephson coupling the shift is still too 
small to account for the difference in the transition tem- 
peratures between Bi-2212 bulk und thin film samples. 

Finally we have shown that the renormalization group 
approach fails for intermediate Josephson coupling. The 
reason for this failure has been argued to lie in the one- 
dimensional nature of the Josephson vortex strings. One 
would welcome a renormalization group theory for this 
system, but if and how it may be constructed remains a 
question for the future. 
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APPENDIX A: 



As discussed in Sec. |lll| we have to integrate the expression exp {—(3k |r< — Iq| — (3 k \r\ — r 3 -|) over and r,-. The 
integral over is given by 

Ji = / d 2 n exp (-(3k |r< -ri\-0KT- (3k \y{ -r 5 \+ (3k \v{ - rj,|) 

Jr^lri-rj \<T+dr ^ ' 



cos 9 



= rdrexp (^—(3kt — (3k \rj — r{| + (3k |r{ — Iq|^ 
x J d9exp^-(3K \rj - + 2— — 

Expansion for small r/\Tj — r z Q \ yields 

Xi = t dT exp ^ — (3kt — (3k \tj — r^| + (3k |r^ — r ol^ J ^0 exp 

= 2nTdTexp (-(3k \vj - t % q\ - (3kt - (3k \vj - rj| + (3k |rj - r l \j I q {(3k), 
and integration over Vj finally gives 

/ d 2 r ] l 1 = 2nrdT exp (-0kt + (3k |rj - r*|) I {(3k) \r{ - rj,| 2 K 2 ((3k \t{ - rj|) 



^ki-rjl 1 



(Al) 



cos 8 



(A2) 



(A3) 



This is the result stated above. 
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FIG. 1. The trajectories of Eqs. fl2.24|) and ( |2.28; ), obtained 
by numerical integration for z(0) = 0.03 and several values of 
K(0). The flow is to the left. 



FIG. 2. The three distinct configuration of a small vor- 
tex-antivortex pair with respect to the JV's. The arrows de- 
note pancake vortices and the thick gray lines JV's. See text. 



FIG. 3. The trajectories for a system with strong Joseph- 
son coupling, but vanishing electromagnetic coupling between 
the layers. The flow is to the left. The dashed curve would 
be the critical trajectory in the original BKT case. 



FIG. 4. Definition of vortex coordinates. The circles de- 
note pancake vortices and the dashed line is a JV. The arrows 
point in the direction of the magnetic field within the JV. 



15 



